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Polarized proton beams at HERA can currently only be produced by extracting a beam 
from a polarized source and then accelerating it in the three synchrotrons at DESY. In 
this paper, the processes which can depolarize a proton beam in circular accelerators 
are explained, devices which could avoid this depolarization in the DESY accelerator 
chain are described, and specific problems which become important at the high energies 
of HERA are mentioned. At HERA's high energies, spin motion cannot be accurately 
described with the isolated resonance model which has been successfully used for lower 
energy rings. To illustrate the principles of more accurate simulations, the invariant 
spin field is introduced to describe the equilibrium polarization state of a beam and the 
changes during acceleration. It will be shown how linearized spin motion leads to a 
computationally quick approximation for the invariant spin field and how to amend this 
with more time consuming but accurate non-perturbative computations. Analysis with 
these techniques has allowed us to establish optimal Siberian Snake schemes for HERA. 

1. INTRODUCTION 

In contrast to polarized high energy electron beams which can become polarized by 
the emission of spin flip synchrotron radiation, proton beams do not become polarized 
after acceleration. High energy polarized protons can only be produced by accelerating 
a beam from a polarized ion source. The highest energy reached so far was 25GeV in 
the AGS and a lot of care had to be taken during acceleration to preserve this injected 
polarization [|JBJ. To explain depolarization in circular accelerators, some concepts from 
spin dynamics have to be introduced. When a particle with charge q moves through a 
magnetic field, the motion of the classical spin vector in the instantaneous rest frame is 
described by the Lorentz force equation and the Thomas-BMT equation 

q -{B ± }xp, ^ = -J-{(Gi + l)B ± + (G + l)B«}x3, (1) 



dt rwy dt rwy 

where B± and B± are the magnetic field components perpendicular and parallel to the 
particle's momentum p . At high energy where (G7 + ~ G, the spin motion is 
independent of energy and a fixed field integral of 5.48Tm leads to a spin rotation of 
180° , in contrast to the orbit deflection which varies with I/7. For fixed orbit deflections 
and thus fixed ratio of B±/j, the spin precession rate increases with energy as we now 
describe. 
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In purely transverse magnetic fields, the Thomas-BMT equation has the same structure 
as the Lorentz force equation up to a factor G7 + 1 . The spin therefore rotates like the 
momentum but with a magnified rate. At the HERA energy of 920GeV, the magnification 
factor is G7 = 1756 so that the spin is rotated by 100° when a proton's direction is 
altered by lmrad in a transverse magnetic field and the spin rotates 1756 times around the 
vertical while a particle makes one turn on the design orbit of a flat circular accelerator, 
where the fields are vertical. The number of spin rotations which a particle performs while 
it travels along the closed orbit once is called the spin tune u . 

When the spin tune is integer, the spin comes back to a field imperfection with the same 
direction after one turn and the effect of the field error can add up coherently from turn 
to turn. This resonant depolarization at integer spin tunes uq is called an imperfection 
resonance @. 

When viewed at a fixed azimuth 9 of the accelerator, the particles appear to perform 
harmonic oscillations around the closed orbit with the frequencies v x , v y , and v T for 
horizontal, vertical, and longitudinal motion. These are called the orbital tunes. Some 
of the fields through which a particle propagates will therefore oscillate with the orbital 
tunes. Whenever the non-integer part of the spin tune is equal to plus or minus one of 
these frequencies, the resulting coherent perturbation can lead to depolarization. The 
coherent depolarization at the first order resonance condition Vq — m + v%. is called an 
intrinsic resonance 0. Here the notation Ux—u x , Vi=-v x , v%=v yy va=-v Vi p^,=v t , v§=-v t 
is used. Since the spin tune changes with energy (in a flat ring v§ = G7) resonances will 
have to be crossed at some energies during acceleration. 

After one turn around the accelerator, all spins of particles on the closed orbit have 
been rotated by 2ttu around a unit rotation vector n . This vector is determined by 
the accelerator's main guide fields and small field imperfections only perturb Hq weakly 
except at energies where the guide fields would produce an integer spin tune. Then, when 
viewed from a fixed azimuth, spins would come back after one turn apparently without 
a rotation. Close to imperfection resonances the remaining rotation and therefore the 
direction of n will be dominated by the influence of field errors. While the energy changes 
during acceleration, no changes its direction strongly at these resonances. Whenever this 
change is sufficiently slow, spins which are initially parallel to no will follow the change of 
f?o adiabatically. Imperfection resonances can therefore be crossed either by making the 
field imperfections small enough or by making them so strong that n already starts to 
get influenced by the field errors sufficiently long before the resonance and then changes 
slowly enough to let all spins follow adiabatically while the resonance is crossed. Special 
magnets for enhancing this effect without disturbing the orbit are referred to as partial 
snakes. So far solenoid magnets have been used Q but for the AGS a helical dipole partial 
snake is under construction JJ. 

The motion of spins along phase space trajectories is dominated by the main guide fields 
on the closed orbit except close to an intrinsic resonance, where the coherent perturbations 
described above can dominate over the main guide fields. When the emittance of the beam 
and therefore the amplitude of the perturbations is sufficiently small, intrinsic resonances 
can be crossed without loss of polarization. Polarization in the core of the beam will 
therefore be only weakly influenced when crossing intrinsic resonances. If a strong coherent 
perturbation is slowly switched on and off, an effect similar to adiabatically following n 
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occurs and polarization is conserved. Therefore, while an intrinsic resonance is crossed, 
perturbations influencing particles in the tails of a beam will slowly increase already before 
the resonance and this adiabatic conservation of polarization can occur. In intermediate 
parts of the beam, however, the polarization is lost. This type of depolarization can be 
overcome by slowly exciting the whole beam coherently at a frequency close to the orbital 
tune which causes the perturbation. All spins then follow the adiabatic change of the 
polarization direction and the resonance can be crossed with little loss of polarization. 
The excitation amplitude is then reduced slowly so that the beam emittance does not 
change noticeably during the whole process. This mechanism has recently been tested 
successfully at the AGS [0. An older technique of avoiding depolarization at strong 
intrinsic resonances utilizes pulsed quadrupoles to move the orbital tune within a few 
microseconds just before a resonance so that the resonance is crossed so quickly that the 
spin motion is hardly disturbed ||. 

So far no polarized beam has been accelerated to more than 25GeV ||. But the pos- 
sibility of polarized proton acceleration has been analyzed for RHIC (250GeV), for the 
TEVATRON (900GeV), and for HERA (920GeV). When accelerating through approx- 
imately 5000 resonances in the case of HERA, even very small depolarization in every 
resonance crossing would add up to a detrimental effect. 

It was mentioned below equation ([!]) that in a fixed transverse magnetic field the deflec- 
tion angle of high energy particles depends on energy, whereas the spin rotation does not 
depend on energy. It is therefore possible to devise a fixed field magnetic device which ro- 
tates spins by tt whenever a high energy particle travels through it at the different energies 
of an acceleration cycle. Such field arrangements which rotate spins by it while perturbing 
the orbit only moderately are called Siberian Snakes [[J. The rotation axis is called the 
snake axis and the angle of this axis to the beam direction is referred to as the snake angle 
ip . Let us consider a Siberian Snake with snake angle ipi at one point in a flat ring and 
a second Siberian Snake with snake angle ?/> 2 at the opposite side of the ring where the 
spin has rotated by G^y/2. The spin rotation around the vertical between the Siberian 
Snakes is described with Pauli matrices by the quaternion cos(7rG7/2) + i sin(iiG'y /2)a 2 . 
The rotation by the first Siberian Snake is described by i[sm(i/;i)ai + cos(^i)o"3] . The 
total rotation for one turn around the ring is then described by 

i [sin (■?/>! )<7i + cos(^i)cr 3 ] • [cos(7rG7/2) + i sin(7rG7/2)<7 2 ] 

• i[sin(^ 2 )c r i + cos(^ 2 )o"3] • [cos(7rG7/2) + isin^G-y /2)a 2 ] 
= i[sin(V>i + 7rG7/2)(Ti + cos(V>i + nG-f/2)a 3 ] 

■ i[sm(i/j 2 + irG'y/2)cr 1 + cos(-?/> 2 + nG'y /2)cr 3 ] 
= -cos(V>i — 1P2) +isin(ipi — 1^2)^2 ■ (2) 

For ipi — ip2 = 7r/2 the spins rotate in total 1/2 times around the vertical no during 
a complete turn around the ring, giving vq = 1/2. All imperfection resonances and, 
since the orbital tunes cannot be 1/2, also all first order intrinsic resonances are avoided 
by the insertion of these two Siberian Snakes, and polarized beam acceleration to very 
high energy could become possible. Siberian Snakes can only be used at sufficiently 
high energies since their fields are not changed during acceleration of the beam and they 
produce orbit distortions which are too big for energies below approximately 8GeV |i0| . 
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2. THE DESY ACCELERATOR CHAIN FOR POLARIZED PROTONS 

For HERA a polarized proton beam would be produced by a polarized H source. Then 
it would be accelerated to 750keV in an RFQ and then to 50MeV in the LINAC III from 
where it would be accelerated in the synchrotron DESY III to 7.5GeV/c. In the next 
ring, PETRA, 40GeV/c are reached, and HERA finally accelerates to 920GeV/c. The 
four main challenges for obtaining highly polarized beams in HERA are: (1) Production 
of a 20mA pulsed H beam. (2) Polarimetry at various stages in the acceleration chain. 
(3) Acceleration through the complete accelerator chain with little depolarization. (4) 
Storage of a polarized beam at the top energy over many hours with little depolarization. 

Polarized protons are produced either by a polarized atomic beam source (ABS), where 
a pulsed beam with 87% polarization for 1mA beam current has been achieved, or by an 
optically pumped polarized ion source (OPPIS), where pulsed beams with 60% for 5mA 
have been achieved. Experts claim that 80% polarization and 20mA could be achievable 
with the second type of source. The current source at DESY produces 60mA but the 
maximal current of 205mA in DESY III can already be achieved with a 20mA source. 

Polarimeters will have to be installed at several crucial places in the accelerator chain. 
The source would contain a Lyman-a polarimeter |i~Tfl . Another polarimeter could be 



installed after the RFQ |i~2| . This could not be operated continuously since it disturbs 



the beam. The transfer of polarized particles through the LINAC III could be optimized 
with a polarimeter similar to that in the AGS LINAC; and like the AGS, DESY III could 
contain an internal polarimeter ||. Polarization at DESY III energies has been achieved 
and measured at several labs already. It is different with PETRA and HERA energies; 
for these high energies there is no established polarimeter. Here one has to wait and see 



how the novel techniques envisaged and developed for RHIC will work flT"3 |. 

Since DESY III has a super period of eight, only 4 strong intrinsic first order resonances 
have to be crossed. They are at values for the spin tune G7 of 8 — u y , + v y , 16 — u yi and 
8 + v y . Depolarization can be avoided by jumping the tune with pulsed quadrupoles in a 
few microseconds or by excitation of a resonance with an RF dipole. A solenoid partial 
snake would be used to cross the one strong imperfection resonance at G7 = 8. All these 
methods have been tested successfully at the AGS and it is likely that a highly polarized 
proton beam could be extracted from the DESY III synchrotron at 7.5GeV/c . 

In PETRA it would be very cumbersome to cross all resonances which can be seen 
in figure [l] (middle). Since Siberian Snakes can be constructed for the injection energy 
of PETRA [|14]] it will be best to avoid all first order resonances with two such devices. 
There is space for Siberian Snakes in the east and the west section of PETRA. 

3. SPECIFIC PROBLEMS FOR THE HERA RING 

HERA is a very complex accelerator and a brief look already indicates four reasons why 
producing a polarized beam in HERA is more difficult than in an ideal ring. (1) HERA has 
a super periodicity of one and only an approximate mirror symmetry between the North 
and South halves of the ring. Therefore more resonances appear than in a ring with some 
higher super periodicity and special schemes for canceling resonances in symmetric lattices 



[|nj cannot be used in such a ring. (2) The proton ring of HERA is on top of the electron 



ring in the arcs, and the proton beam is bent down to the level of the electron ring on 
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both sides of the three experiments HI, HERMES, and ZEUS in the North, East, and 
South straight sections. The HERA proton accelerator is therefore not a flat ring. The 
destructive effect of the vertical bends can, however, be eliminated by so called flattening 
snakes p6|JT7| which let the spin motion in pairs of vertical bends cancel and makes 



n vertical outside the non-flat sections of HERA. (3) There is space for spin rotators 
which make the polarization parallel to the beam direction inside the collider experiments 
while keeping it vertical in the arcs, and there is also space for four Siberian Snakes. 
But installing more than four Siberian Snakes would involve a lot of costly construction 
work. Simulations have shown that 8 snakes with properly chosen snake angles would be 
desirable. However, if one does not choose optimal snake angles, then four-snake-schemes 
can be better than eight snake schemes ||18|| . (4) The energy is very high and therefore 
the spin rotates rapidly. If HERA had been designed for polarized proton acceleration, 
several parts of the ring would probably have been constructed differently. 



4. APPLICABLE THEORY AND SIMULATION TOOLS 

4.1. The isolated resonance model 

In the isolated resonance model, the field components which perturb the spin of a 
particle that oscillates around the closed orbit are Fourier expanded. The perturbation of 
spin motion is then approximated by dropping all except one of the Fourier components. 
When z describes the phase space coordinates relative to the closed orbit and 9 describes 
the accelerator's azimuth, the Thomas-BMT equation ([]]) has the form dsj d9 = Q(z, 9) x s. 
The precession vector Q can be written as f2o(#) + dj{z,9) with a part on the closed 
orbit and a part which is linear in the phase space coordinates z. For spins parallel to 
the rotation vector on the closed orbit n*o{9) only the components of lj(z,9) which are 
perpendicular to n*o perturb the polarization. We now choose two mutually orthogonal 
unit vectors m and Zo which are perpendicular to n and precess around flo according to 
the Thomas-BMT equation on the closed orbit. The frequency of their rotation is given 
by the spin tune z/ . 

In this model a depolarizing resonance occures when a Fourier component of uj(z(9), 9) 
rotates with the same frequency as rho and Iq so that there is a coherent perturbation of 
the spins away from no- The Fourier component e VQ for this frequency is obtained from 
the Fourier integral along a trajectory z(9), 

e vo = Hm — - / u{z{9), 9) ■ (m + il Q )dO . (3) 

These resonance strengths are shown in the figure [l] (top), (middle), and (bottom) for the 
three proton synchrotrons at DESY. They were all computed for an oscillation amplitude 
of z(9) corresponding to the one sigma vertical emittance of 47rmm mrad. 

4.2. The invariant spin field 

Already at extraction from PETRA the polarized beam would have somewhat more 
energy than any other polarized proton beam so far obtained and one has to ask whether 
the isolated resonance model successfully used so far for describing depolarization is still 
applicable. To understand whether the isolated resonance model describes spin motion 
at HERA accurately, we introduce the invariant spin field of a circular accelerator. It has 
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been mentioned that a particle on the closed orbit has to be polarized parallel to no in 
order to have the same polarization after every turn. Similarly, one can ask if the whole 
field of spin directions for particles at different phase space points can be invariant from 
turn to turn. 

Each particle can have a different spin direction at its phase space point z and each 
of these spins propagates with a different precession vector Q(z(9),6) in the Thomas- 
BMT equation. A spin field n(z) which is invariant after one turn around the ring is 
called an invariant spin field or a Derbenev-Kontratenko n-axis ||. A beam which is 
polarized parallel to this invariant spin field at every phase space point does not change 
its polarization state from turn to turn. Particles change their location in phase space 
from some initial phase space coordinate z% in the Poincare section at azimuth 9 to some 
final coordinate after one turn Zf = M(zj) according to the one turn map. And spins 
change their directions according to the one turn spin transport matrix R(zi), but the 
invariant field of spin directions n(zi) does not change after one turn. This requirement 
is encompassed by the periodicity condition 

n{M{z i ))=R{z i )n{z i ) . (4) 



Note that the polarization state of a particle beam is in general not invariant from turn to 
turn when all particles are initially completely polarized parallel to each other, but rather 
when each particle is polarized parallel to n(z) at its phase space point z. In this case 
the polarization of a particle will be parallel to n(zi) whenever it comes close to its initial 
phase space point Z{ during later turns around the ring, as long as n{z) is sufficiently 
continuous. When two particles travel along the same trajectory, the angle between their 
two spins does not change. When a particle is initially polarized with an angle <fi to n(z), 
it will therefore be rotated around n(z) every time it comes close to Zi, but it will still have 
the angle to the invariant spin field. The time averaged polarization at ^ will therefore 
be parallel to n(zi), but it can only have the magnitude 1 if the spin was initially parallel 
to the invariant spin field. However, even if all particles are initially polarized parallel to 
n(z) , the beam polarization is not 1 but < n > where < . . . > denotes an average over the 
beam. The maximum average polarization that can be stored in an accelerator at a given 
fixed energy is therefore | < n > |. It was first pointed out in |19| that this maximum 
polarization can be small in HERA. 

Since the spin dynamics depends on energy, the invariant spin field n(z) will change 
during the acceleration process. If this change is slow enough, spins which are parallel 
to n(z) will follow adiabatically. However, if the change is too rapid, polarization will be 
lost. It is therefore good to have < n > close to 1 not only at the collider energy but 
during the complete acceleration cycle. Four problems occur when the different directions 
of n(z) are not close to parallel for all particles in the beam. (1) Sudden changes of n(z) 
reduces the polarization. (2) The average polarization available to the collider experiment 
is reduced. (3) The polarization involved in each collision process depends on the phase 
space position of the interacting particles. (4) Measuring the polarization in the tail of 
the beam will not give accurate information on the average polarization of the beam. 
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4.2.1. Linearized spin orbit motion 

At azimuth 9, a spin can be described by a usually small complex coordinate a with 



s — ^{a}m (9)+^s{a}l (8) + ^l — \a\ 2 n (8). When the spin coordinates a and the phase 
space coordinates are linearized, one approximates an initial spin by Si ~ 9?{a;}m (0) + 
S{a!i}Zo(0) + no(0) at azimuth and the final spin after one turn around the accelerator 
by Sf = 5J{a/}m (0) + $5{a/}Z (0) + n (0) where «/ is determined by the 7x7 one turn 
transport matrix M_ 77 , 



(5) 



whereby M is the 6x6 dimensional one turn transport matrix for the phase space variables, 
the exponential describes the rotation of the spin components a by the spin tune z/ around 
n , the row vector G T describes the dependence of spin motion on phase space motion 
to first order, and the 6 dimensional zero vector shows that the effect of Stern Gerlach 
forces on the orbit motion is not considered. 

We now write the components perpendicular to n of the invariant spin field as a com- 
plex function n a (z) and use a 7 dimensional vector fti to obtain the first order expansion 
of n(z). The linearized periodicity condition for the invariant spin field is 



n 1 (z)=[ , n^Mz) = M^z) . (6) 




This equation can be solved for fti after the matrices are diagonalized. Let AT 1 be the 
column matrix of eigenvectors of the one turn matrix M_. The diagonalized matrix of 
orbit motion A = A M_ A -1 has the diagonal elements exp(z27rz/fc) given by the orbital 
tunes vi = u x , u 2 = —v x -, etc. We now need the 7x6 dimensional matrix T which is the 
column matrix of the first 6 eigenvectors of M_ 77 and has the form 




T= [- Br , TA = M 77 T, (7) 



where the 7th components of the eigenvectors form a vector B. If a linear function fii(z) = 
Kz of the phase space coordinates can be found, which satisfies the periodicity condition 
(|6|), then an invariant spin field has been determined. Inserting the form z\ = Kz into 
equation (|6|) and multiplying the resulting condition K_ M_ = M 77 K by A' 1 from the 
right leads to K A~ X A = M 77 KA ~ l . Therefore K A' 1 is the 7 x 6 dimensional matrix of 
eigenvectors T satisfying equation (|7|) and we conclude that there exists a unique linear 
invariant spin field given by 

ni(z) = TAz. (8) 

In the linear approximation of spin motion, the invariant spin field is simply computed 
via the eigenvectors of the 7x7 spin orbit transport matrix. This matrix M 77 can 
be computed in various ways, for example by multiplying the individual spin transport 



matrixes of all elements [20] or by concatenating spin transport quaternions of individual 



elements as done in the program SPRINT In the normal form space belonging to 



8 



the diagonal matrix A the coordinates are given by the actions Jj and the angle variables 
$j with 

Az = . (9) 

The average over all angle variables of a phase space torus then leads to the average 
opening angle of 



< 0(n, n ) >^ atan(^/< \n a \ 2 >) = atan(, ^2(\B 2k -i\ 2 + |-B 2fc | 2 )J fe ) , (10) 

\ k=l 

where the are the 7th components of the eigenvectors in equation (|7|). 

These opening angles are shown for DESY III in figure |2] (top) and it is apparent that at 
the places where resonant spin perturbations are described by a large resonance strength, 
the invariant spin field has a large opening angle. It is obvious when comparing with 
the resonance strength of figure [l] (top) that the influence of different resonances does 
not overlap in the linearized spin approximation. At PETRA energies of up to 40GeV, 
the resonances already come very close to each other as seen when comparing figure 
|2] (middle) with figure [1] (middle) and one can only barely expect an isolated resonance 
approximation to lead to accurate results. For high energies between 780 and 820GeV/c 
in HERA, figure [] (bottom) clearly shows that one cannot speak of isolated resonances. 
Often the influences of 4 resonances overlap. 

The approximation of linearized spin motion contains all first order orbital frequencies, 
since it linearized the precession vector Q with respect to z. However, in contrast to the 
isolated resonance model, none of these resonances is ignored and the effect of overlapping 
resonances can be seen. 

It is possible to recover the first order isolated resonance strength from the one turn 
spin orbit transport matrix. In analogy to the complex notation for the spin component 
perpendicular to n*o, the perturbing precession vector uj is expressed by a complex function 
uj(z, 9) as uj = ^ft{uj}rho + Q{uj}Iq + (ul • n*o)n*o- Inserting this into the Thomas-BMT 
equation, one obtains 



a' = —i\J\ — a 2 uj + ia{uj ■ no) • (11) 

In the case of spins which are nearly parallel to no, one can linearize in a and z. For a 
spin which was initially parallel to no one obtains a{9) m — i Jq ujdO . Comparing with 
equation (|) one can express the resonance strength as e UQ = ilimN^^ ^j^a(2iiN) . The 
resonance strength can therefore be computed from Mj 7 /N for large N. The computation 
becomes very efficient if one uses M77 = (M^) 2 iteratively. 

The coordinate vectors m (27r) and Iq(2tt) to which a(27r) refers have rotated by the spin 
tune Uo, whereas the final spin coordinate «/ computed by M 77 refers to the coordinate 
vectors m (0) and Iq(0). Therefore a(27rAQ = a/ exp(— i27iNi> ), and e UQ can be computed 
from powers of the one turn matrix, which can most efficiently be evaluated in diagonal 
form, 

e V0 = 1 lim -^—a(2nN)=i lim -^-(0, e" iAf2 ™°) ( if ., Q ) ( fl (12) 
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i N-1 

= Uim e -w*"*_!_ ^[eW-i-V^&A-WW (13) 

N->oo 2nN n 

3=0 

i N-1 

= te-^G^AkmZm lim — £ e «*iC«*-«*> (14) 

3=0 

where one has to sum over equal indices k, I, and m. This formula shows that the 
resonance strength is always zero, except at a resonance condition v = m + vy.. At such 
a spin tune, the resonance strength is given by 



271-1^=^1 = |G T A~ 1 diag(0...1...0)A^ = \G T A' 1 (0...,Jj k e i ' s ' k ...0) T \ = \G ■ v k \^J k ■ (15) 

The 1 in the diagonal matrix is in position k. Here A~ 1 (0...yQ k ~e l ' s ' k ...0) T is the initial 
value for a phase space trajectory which has only Fourier components with frequencies v k 
plus integers and the fcth eigenvector vj, of M has been used. The infinite Fourier integral 
in equation (^) has been reduced to the scalar product between the bottom vector of M 77 
and an eigenvector of M. This very simple formula is used in the program SPRINT. 

4.3. Non-perturbative methods 

While one does not drop Fourier coefficients in the approximation of linearized spin 
motion, there are other limitations. The approximation is no longer justified when \n a \ 
becomes large, which happens close to resonances in the figures ^. Therefore the validity 
of linearized spin motion had to be be checked by computing the invariant spin field non- 
perturbatively. In the last few years two iterative higher order and three non-perturbative 
methods of computing the invariant spin field have been developed . All of these meth- 
ods agree within their ranges of mutual applicability. The invariant spin field obtained 
from a non-perturbative method contains the effect of all Fourier coefficients in Q. When 
comparing this spin field with n*i, it was found that linearized orbit motion describes the 
opening angle and thus the maximum storeable polarization well in domains where the 
opening angle is small. At the critical energies, where the maximum polarization is low 
during the acceleration process, non-perturbative methods become essential for simulation 
and results obtained with the computationally quick linearization of spin motion should 
always be checked with more time consuming non-perturbative methods if possible. 

One application of this strategy is the filtering method ||23|| . Four or eight Siberian 



Snakes are inserted into HERA to fix the spin tune to 1/2 for all energies and to let no be 
vertical in the flat arcs. These conditions do not fix all snake angles. Currently there is, 
however, no established formula to determine good snake angles. Since the opening angle 
of n(z) is such a critical quantity for high energy polarized proton acceleration, we have 
decided to maximize < n > by choosing snake angles. A computer code was written which 
tested approximately 10 6 snake schemes and computational speed was therefore essential. 
Linearized spin motion was used to find the 8-snake-schemes with smallest average value 
of \n a \ over the acceleration cycle. These filtered snake schemes then also had relatively 
small opening angles when computed non-perturbatively with stroboscopic averaging [^T] . 



Two other indications showed that this filtering leads to good snake schemes. (1) Track- 
ing simulations of the complete ramp process showed that the snake schemes found by 
filtering leads to less depolarization |24| than other schemes which were initially proposed. 
(2) Computation of the amplitude dependent spin tune, which can only be performed when 
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n(z) has been found non-perturbatively, shows that snakes schemes found by filtering have 
significantly less spin tune spread over orbital amplitudes than other proposed schemes 

p5fl . With the optimal scheme for four Siberian Snakes in HERA it turned out to be 
possible to accelerate in computer simulations approximately 65% of the beam to high 
energy with little loss of polarization as long as no closed orbit distortions were present 

EH]. In simulations, the current 1mm rms closed orbit distortions lead to depolarization 



2(J • Therefore either the closed orbit will have to be controlled more accurately or tech- 
niques which make the spin motion less sensitive to closed orbit distortion [^7j will have 
to be utilized. 
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Figure 1: Resonance strength for DESY III, PETRA, and Figure 2: Opening angle for DESY III, PETRA, and HERA. 
HERA. 



